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We study the absorption spectrum for a strongly degen- 
erate Fermi gas confined in a harmonic trap. The spectrum 
is calculated using both the exact summation and also the 
Thomas-Fermi (TF) approximation. In the latter case, rela- 
tively simple analytical expressions are obtained for the ab- 
sorption lineshape at large number of trapped atoms. At zero 
temperature, the approximated lineshape is characterized by 
a f 1 — z 2 ) 5//2 dependence which agrees well with the exact nu- 
merical calculations. At non-zero temperature, the spectrum 
becomes broader, although remains non-Gaussian as long as 
the fermion gas is degenerate. The changes in the trap fre- 
quency for an electronically excited atom can introduce an 
additional line broadening. 
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In recent years there has been a great deal of interest in 
the dilute gas of trapped atoms cooled to temperatures 
below 1 /iK. At such low temperatures, an important 
role is played by the quantum statistics of atoms. Bosons 
tend to occupy the lowest translational level of the trap 
to form the Bose - Einstein condensate On the 

other hand, fermions can still occupy many trap levels 
(predominantly up to the Fermi level) , as it was observed 
recently Q . The physical properties of a Fermi gas (such 
as the specific heat) depend on the number of atoms in 
the system |||| , as well as on the trap anisotropy || . 

The optical spectroscopy has proven to be useful in 
getting information about cold atomic gases Jl|-§,0^O). 
The Fermi-Dirac (FD) statistics is known to change the 
optical response of the system compared to the classical 
one. The signatures of quantum degeneracy emerge in 
the scattering spectra of homogeneous [jl0| and trapped 
|jj Fermi gases. Effects of quantum statistics are also 
featured in the scattering of short laser pulses from a 
trapped Fermi gas Furthermore, the spontaneous 
emission appears to be inhibited in a cold Fermi gas [jll) . 

The effects of quantum degeneracy should manifest in 
the absorption spectra as well. The aim of the present 
paper is to investigate absorption by a cold Fermi gas 
confined in a harmonic trap. The analysis concentrates 
on the degenerate Fermi gas (i.e. very low temperatures), 
for which the quantum statistics of the atoms plays an 
important role. The theory involves exact calculations, 
as well as the Thomas - Fermi (TF) approximation. Con- 
sequently, the absorption spectrum is analyzed both for 
small and large numbers of trapped atoms. 

Consider a gas of Fermi atoms confined in a harmonic 



trap. The harmonic approximation is relevant for the 
traps used in recent experiments j|. We shall neglect 
atomic collisions, since the s-wave collisions are forbidden 
between the spin - polarized fermions. Consequently, one 
can make use of the following one - atom Hamiltonian: 



Hl-at = \g) H g (g\ + ^ \ ex ) (TWO, ex + H e , 



(ex 



(1) 



where \g) and \ex) represent the ground and an excited 
electronic state of an atom, Tiu>Q :ex is the excitation en- 
ergy. It is noteworthy that the ground electronic level 
of the fermion atom has a number of magnetic sublevels 
over which the summation is to be carried out in the 
Hamiltonian ([l]). However, such a summation is not nec- 
essary if the atoms are spin - polarized, as it is the case 
in the experiment by DeMarco and Jin on 40 K atoms [[| . 
Here also H g (H ex ) is the Hamiltonian for the transla- 
tional motion of a trapped atom in the ground (excited) 
electronic state: 



H g<ex (r, p) 
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where p = —iKV is the momentum operator, M is the 
atomic mass, £l g (Sl ex ) is the frequency of the transla- 
tional motion along the x axis if the atom is in the ground 
(excited) electronic state, and the dimensionless param- 
eters A y and A z describe the extent of anisotropy of the 
trap. Note that the frequency f2 g can be generally dif- 
ferent from n ex due to the changes in the magnetic mo- 
ment of the atom following its transition to an excited 
electronic state. The effects related to this fact will be 
explored using the TF approximation. 

The lineshape of the absorption spectrum is given by 



I(u) = J2p*\(f\V\i)\ 2 6(Lu-tu fl ). 



(3) 



are the ini- 



Here \i) = \g)\n) g and |/) = \ex}\m}^ 
tial and final states of an atom, |n) fl = \n x ,n y ,n z ) g 
and |m) = \m Xl m yi m z ) are the atomic translational 
states characterized by the energies 



e 3 X x = TiClg t ex (Tlx + ^ y n y + X z n z ) + eg' 



g,ex 



(4) 



l q,ex 

where £q 



(1 + A y + A z ) /2, and cj fl 
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(e^ — e^) /h is the transition frequency. Here also pi = 
p{ £ a i ) = [exp(/?££ - /3/ji) + lp 1 is the FD distribution 
function for the trapped atoms, fj, is the chemical poten- 
tial and (3 = l/ksT. The operator 



1 



V 



{\ex) d ex exp (inx) (g\ + h.c.} 
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describes the interaction between an atom and the elec- 
tromagnetic field propagating along the x axis, k is the 
wave number of the light, and d ex is the atomic transition 
dipolc moment along the polarization of the light. 

Consider first the absorption spectrum using the exact 
summation over the translational levels. At this stage, 
it is assumed that il g = Q ex = f2, yet the trap can still 
be anisotropic. At zero temperature, only the levels with 
n x + X y n y + X z n z < np = Ep/hQ, are occupied by the 
atoms, where Ep = /i| T=0 is the Fermi energy. In such 
a situation, the absorption lineshape takes the form: 

tip oo 

I ( UJ ) = E \d ex | 2 ^2 ^2 K x 5 (uj - uj , ex - mil) 
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with j x = min (n x , n x + m), where a = k(Ti/2M£1) 1 , 

and the factor K. = Y^^'^ 1 

represents a number of occupied translational states 
\n x , n y , n z ) for a fixed value of n x , the brackets [...] label- 
ing the integer part of a number. The Fermi number np 
is determined by the condition J2n F =o ^ x = where N 
is the number of trapped atoms. The factor K x reflects 
the trap geometry. For traps with a cylindric symme- 
try (X y = X Z = A), one finds K x = ([q x ] + 1) ([q x ] + 2) /2, 
where q x = (np — n x ) / A. For an anisotropic trap of a 
cigar shape (X y ,X z » 1), one has K x = 1 provided the 
number of trapped atoms is small enough (np < X y ,X z ). 
In such a situation, the trap becomes one-dimensional 
(ID), giving np = N — 1. On the other hand, for 
an isotropic three-dimensional (3D) trap (X y = X z = 1), 
one arrives at K x — (np — n x + 1) (np — n x + 2) /2. If 
the number of atoms is sufficiently large (np > X y ,X z ), 
the anisotropic traps of cigar shape (X y , X z 3> 1) are no 
longer one dimensional, since the Fermi energy is then 
greater than the energy of the translational quanta in 
the y and z directions. Such a situation corresponds to 
the recent experiment by DeMarco and Jin 0). 

Figure |l| shows the absorption lineshapes for vari- 
ous degrees of the trap anisotropy in the case where 
X y = X z = X. A single excited electronic state \ex) has 
been taken into account in these and the subsequent fig- 
ures. Oscillations are seen clearly in the thick solid line 
representing a purely ID case (A = 20), as well as in 
the thin one corresponding to an anisotropic 3-D case 
(A = 5). This can be related to the oscillations of the 
density of Fermi atoms in the one-dimensional Jl2 ] and 
anisotropic three dimensional traps Q at a sufficiently 
small number of the trapped atoms. Oscillations do not 
appear in the lineshape of an isotropic 3D trap (A = 1). 



Note that in contrast to a single trapped atom |13|], the 
zero-temperature lineshape of the trapped Fermi gas has 
a cut off at the frequencies smaller than ujo >ex . This can 
be explained by the fact that fermions occupy excited 
translational levels of the trap at T=0 (up to the Fermi 
level), so that optical absorption can be accompanied by 
a decrease in the translational energy of the atoms. 




FIG. 1. Absorption lineshape calculated using Eq.(q) for 
N = 20, T = and a = 8. The singular contributions due to 
the delta functions have been smoothened. 

When N or T is increasing, the behavior of a quantum 
system becomes more similar to that of the classical one. 
To get analytical formulas for the absorption lineshape 
at arbitrarily large values of N and T, we shall make 
use of the semiclassical Thomas-Fermi (TF) approxima- 
tion. In the TF approximation, the state of an atom is 
labeled by the radius-vector r and wave vector k = p/h 
(see e.g. refs. HH). The density of such states in the 
six-dimensional phase space (r, k) is (2tt) 3 . The num- 
ber density of the fermion atoms in the phase space is: 



p(r,k,T) 



1 



(2ir) 3 exp {/3H g (r, ftk) - f3u} 



1 



(7) 



where the chemical potential fi is related to the num- 
ber of trapped atoms via the normalization condition 
/ d 3 rd 3 kp (r, k, T) = N. 

Applying the TF approximation, the lineshape reads: 



/ (w) = \dex I 2 / d 3 rd 3 kp (r, k, T) 5 (. 



(8) 



where w r , k = M (f^ - n 2 g ) (x 2 + Xly 2 + X 2 z z 2 ) /2ft + 
hk x n/M + hn 2 /M + ujQ^ ex is the transition frequency. If 
fig = £l e x, the frequency Lu r .k does not depend on the 
atomic position r, so the lineshape is determined exclu- 
sively by the momentum distribution function: 



n(k,T) = J d 3 rp(r,k,T). 



(9) 
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For T = the distribution function is given by 



n (k, 0) 



87V 
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where K F = (2ME F /hf) is the maximum momen- 
tum of the trapped Fermi atoms at zero temperature, and 

1/3 

Ep = hilg (6X y X z N) 1 is the Fermi energy. Note that 
the momentum distribution is isotropic even though the 
trap is anisotropic Q . This leads to an isotropic absorp- 
tion lineshape in the case where Q ex = Q g . Applying the 

,2VV2 



distribution function (|10|), one arrives at the 
behavior of the lineshape if £L ex = Sl g = Q: 



16iV 
5ttA ^ 

ex 



\d ex \ — Umax) /A 2 



1 



5/2 



(11) 



where the central frequency tj max = wq. e 



2 fl is shifted 



by the recoil frequency uo rec — a fi as compared to wo.e 



A = aQ (6NX y X z ) 



1/6 



(12) 



being the spectral half-width. In the experiment by De- 
Marco and Jin [Q using trapped 40 K atoms, f2 = 2ir x 19 
Hz, N = 7xl0 5 and A w 7, giving a rj 36 for w = 4xl0 16 
Hz. Consequently one has A w 6 x 10 5 Hz. This is less 
than the typical radiative linewidths for free atoms. Yet, 
for trap ped fermions the spontaneous emission is sup- 
pressed |ll|, so the Doppler broadening can be dominant. 

The approximated lineshape (|ll ) depends on the trap 
anisotropy exclusively through the characteristic fre- 

1/3 

quency Q c har = n(X y X z ) which is a measure of the 
trap hardness. The bigger VL c har is, the tighter is the 
trap, and the broader is the absorption spectrum. In 
fact, the maximum momentum of the atoms is larger in 
tighter traps (for the same number of trapped atoms) 
leading to the increase in the Doppler broadening. For 
instance, compared to an isotropic trap (X y = X z = 1), 
the spectrum of a squeezed trap (X v , X z > 1) is broader. 

Exact and approximated lineshapes are presented in 
Figs. U and [| For an isotropic trap (Fig. ||), the agree- 
ment appears to be very good, even though the number 
of atoms N is rather small. Deviations are seen only in 
the tails of the spectrum corresponding to the periphery 
of the fermion cloud. In such an area, the fermion density 
becomes small and the TF approximation fails ||. For 
anisotropic traps, the exact spectrum undergoes some os- 
cillations about the approximated one even for relatively 
large values of N (see Fig. ||). In fact, the energy of trans- 
lational quanta depends now on the specific directions of 
atomic motion, so a larger number of trapped atoms is 
needed to populate substantially the translational levels 
in all three directions. 




(co - (fl„ )/£2 

1 0,ex' 

FIG. 2. Absorption lineshape at T — calculated using 
the exact summation (dashed line) and the TF approximation 
(solid line) for a = 9, N = 35 and X x = X y = A = 1. 




50 100 

(co - co„ )/n 

1 0,ex' 

FIG. 3. The same as in Fig. |, but for N = 196 and A = 5. 

Consider next the situation where T ^ and the fre- 
quencies fig and £l ex are not necessarily equal. The line- 
shape (||) takes then the form for an isotropic trap: 

I j |2 poo 

ex r u 

-y 2 -p (w - u) m ax + m ex y 2 /f3h) 2 j , (13) 
where m ex = 1 — f^/^ 2 !, p = (3h/4a 2 fl and f2 = fi 5 . 



If Q g = ilex, the result (13) can be extended readily to 



anisotropic traps. In such a case, the lineshape ( |l3| ) ac- 
quires an extra factor l/X x X y , and the chemical potential 
/i depends on A^Aj,, in addition to T and N. 

We are interested primarily in the strongly degenerate 
Fermi gas ((3Ep ^> 1), for which the Sommerfeld expan- 
sion holds for the chemical potential: 



(14) 



3 



In the opposite (non-degenerate gas) limit ((3Ep — ► 0), 
one has [i — In (f3Ep) 3 /6 , and the lineshape Jl3| ) 



reduces to the Gaussian form if Q a 



kT/E F = 
kT/E F = 0.1 
kT/E f = 0.125 
kT/E r = 0.25 
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FIG. 4. Absorption lineshape at various temperatures for 
an isotropic trap with Q s = Q ex , N = 10667 and a = 9. 



tronically excited atoms leads to an increase and broader 
distribution of frequencies of the optical transitions. 

In summary, we have studied the absorption spec- 
trum by a cold gas of Fermi atoms using both the ex- 
act summation and also the Thomas - Fermi approxima- 
tion. Oscillations have been obtained in the absorption 
spectrum calculated numerically for one-dimensional and 
anisotropic three-dimensional traps at a sufficiently small 
number of trapped atoms and T = 0. No such oscilla- 
tions appear for the isotropic three-dimensional traps. 
Applying the TF approximation, relatively simple ana- 
lytical expressions have been obtained for the lineshape 
of three dimensional traps at a sufficiently large number 
of trapped particles. At zero temperature, the approxi- 
mated spectrum is characterized by a (l — z 2 ) 5 ^ 2 depen- 
dence. At non-zero temperature, the spectrum becomes 
broader, although remains non-Gaussian as long as the 
fcrmion gas is degenerate. The changes in the trap fre- 
quency for an electronically excited atom can introduce 
an additional line broadening. 

The authors would like to express their appreciation for 
the useful discussions with H. Carmichael and T. Moss- 
berg. This work has been completed during the Fulbright 
Scholarship by G. J. at the University of Oregon. 
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FIG. 5. Absorption lineshape for an isotropic trap with 
N = 10667, a = 9, k B T = 0.25£ F and various Q. ex /Q. g . 

Figure |J shows the temperature dependence of the 
lineshapes for Q g = Q ex . At very low temperature 
(kT/Ep = 0.1), the absorption spectrum is seen to be 
close to the zero-temperature limit. As the temperature 
increases, the absorption lineshape becomes broader and 
is no longer characterized by the (l — z 2 ) 5 ^ 2 behavior. 
Yet, the lineshape is still non-Gaussian, since the gas 
is strongly degenerate. Figure ^ shows the absorption 
spectrum for various values of m ex = 1 — fJ^/fi 2 . For 
Q ex > ilg, the maximum position of the spectrum is 
shifted to larger frequencies. For Q, ex < fl g , one has the 
opposite. Furthermore, one can see the obvious increase 
in the spectrum width if fl ex > This is due to the fact 
that an increase in the translational frequency of the elec- 
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